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A*  ALOORITHH  FOR  tFATT-SQUABK  RAT1CWAI  APTWOXIMATIOII  Of  COMPLEX 

Ftt.vnorrn  ci  ti  utiit  circle 


Bo  Jnnra<  n 


Abstract .  An  algorithm  for  rational  approximation  on  the  unit 
circle  will  be  given.  An  irr.<.  li-tti?  application  is  computation  of 
’•ecuraivc  fillin'.  for  use  on  n  digital  computer. 

1.  Definitions 

Let  f(u)  Vc  a  given  ccnplt  .-valued  function  defined  for  real 
values  of  w»  find  let 
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be  a  complex- valued  rational  approximating  function  defined  for  real 

values  of  u>.  The  sets  A  *  {aQ,  a^,  ...  ,  n^}  and  B  *  {b^,  b^,  ...  ,  bQ} 

are  real  parameters  of  R  (oj),  In  the  formulas  below  a.  and  b.  are 

n*n  l  l 

assumed  to  be  equal  to  zero  if  they  do  not  belong  to  A  or  6,  i.e.  if 
their  subscripts  are  outside  the  intervals  0  -  n  or  0  -  n  respectively. 
We  will  always  assume  that  n  >  n. 


2.  Least-squares  approximation 

Wc  wont  to  determine  the  parameters  aQl  a^,  ...  ,  a^  end 
bQ,  bp  ...  ,  bR  such  that 


f(w)  -  R  (w) 
ra,n 


Lg(f  -  R) 


(2) 


attains  n  r.inir.u*.  (.isjunini*  it  cxiutu).  The  "linear"  Approach  of  taking 
the  partial  derivatives  with  ••ccpcct  to  all  the  parameters  is  not 
fcaaibl*  ,  re  it  will  result  in  a  syoten  of  non-linear  equations  which 
is  hard  to  solve.  W-  will  avoid  this  difficulty  by  using  an  algorithm 
giver,  by  the  author  (Jans son,  19<j5). 


3.  Thv  nl/v  r  i  (  h.r.. 

Following  tb.  algorithm  mentioned  in  the  last  paragraph  wo  start 

by  adding  i  real  constant  h  to  f (w),  thus  having  g(u)  ■  f(w)  ♦  h. 

Let  B°  *  {l,  0,  •••  ,0,  .  *.  ,  o}  and  consider  L^( 6  -  R)  as  a  function 

of  aQ,  a^,  ...  ,  hq  only.  We  enn  minimize  Lg(g  -  R)  by  minimizing 
2 

L?(g  -  R)  |  and  this  makes  the  analysis  simpler.  It  may  be  shown  that 
n  unique  minimum  exists.  Furthermore,  it  is  clear  that  Lg(g  -  R)  does 
not  have  a  maximum  value,  for  we  con  place  R(w)  as  for  away  fren  f(x) 

P 

os  we  please.  Since  L^(g  -  R)  is  a  differentiable  function  of 
aQ,  a^,  ...  ,  aQ,  the  minimum  must  occur  vhere  all  of  the  partial 
derivatives 


3  L2U  -  R)2 
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ore  zero. 

From  the  fact  that  h^(c  -  R)  is  a  strictly  convex  function  of 
aQ,  o^,  ...  ,  a^,  it  follows  that  there  is  only  one  set  of  parameters 
for  which  all  of  the  partial  derivatives  are  zero. 

We  have 


I  -  |g(w)  -  \tn(u)|2  " 


■  |g(w)|2  * 
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where 


(fl0 +  *!  +  **»  *ah)>>‘^(nu  al‘fal  *2*  Vl 

(K.  +  ■*■  .. ,  ♦  b^)  ♦  t(bn  b.  ♦  b.  b5+  b 
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bn)  COBU+.  ..♦2(bQ  bfl) 


end 

Be{g(w)  Rn>n(w)}  * 

■  |lta{g(w)f  Rc{Pm(w)  Q^Cu)}  - 

-In  {s(tai)}  In  (Pm(u.)  QOTl]/|Qtl(W)|2. 

Wo  have  that 

pm<">  Q"W  •  aQ  bn  eln“  +  aQ  b^  el(n-1)“  +  ...  ♦ 

+  (*0  Val  V*"  *  Vn)  +  (al  *  —  +  Vlbn>  e'l“ 

♦  ...  +  a  b  e“*mw«  (5) 

m  o 

It  is  now  easy  to  compute  the  partial  derivatives 
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-  Rc{g(w) }  f  b  cos  j  «  ♦  Im  {g(«)}  l  b.  .  sin  j  u 
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cos  n  u 


cos  n  u 


(6) 


Getting  nil  the  part ini  derivative*  equal  to  tero  and  substituting 
into  (3)  vc  get  the  matrix  equation 


M.  •  A°  ■  <], 
A  A 


(7) 


0  . 


is  u  n  x  n  Tooplitz  matrix,  A  is  the  column  vector 
T 

(a  ,  »  •••  •  1  )  •  and  G..  n  column  vector.  Let  us  write  out  the 

0*  1*  *  m  Ai 

elements  of  M.  and  G. 
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To  solve  for  A  wo  just  have  to  compute 


A0  «  M”1  G 
A  A 

which  is  very  easy  since  M  is  a  Toeplitz  matrix  that  can  be  inverted 
.  2 

m  cm  operations.  Here  c  is  a  constant  and  n  is  the  number  of  rows 
and  columns.  An  ulgorithn  for  this  purpose  has  been  given  by 
Robinson  (1963). 

We  have  now  done  the  first  step  and  computed  a  trial  solution  A®, 
Let  us  in  the  next  stop  compute  a  trial  solution  B1.  Wo  will  do  this  by 
minimizing 
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Let  uo  write 
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The  parameters  B  enter  linearly  in  the  expression  above  so  that  ve 
may  take  the  partial  derivatives  and  solve  for  just  as  we  solved 
for  A  . 

Again  ve  get  a  matrix  equation 


*8  fil  *  °B 


(11) 


where  Mg  is  a  Tocplitz  matrix  and  Gg  is  a  column  vector.  We  solve  for 
B^  by  computing 


b1  •  m;1  ob. 


The  elements  of  M  and  G~  may  be  written 
£  & 
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We  hr. vc  now  complete  vl  otu*  iteration  of  the  algorithm  and  we  just  hare 
to  continue  to  dote  mine  A‘  and  B°,  A  and  L  etc  until  we  ore 
close  enough  to  the  point  in  the  parameter  space  AUB  where  the 
minimum  (if  it  exit;.:-)  is  located.  We  have  then  computed  a  sufficiently 
good  approximation  rf  (u>)  „o  6(u).  Obviously 


Hf  (w) 

m,n 


S*U>  -  h  q£(w)  P^(u) 

Q^(w)  Q®(w) 
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is  on  equally  rood  approximation  to  f(u>). 
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